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We uniquely determine the infrared asymptotics of Green functions in Landau 
Or 

gauge Yang-Mills theory. They have to satisfy both, Dyson-Schwinger equations 



and functional renormalisation group equations. Then, consistency fixes the relation 
between the infrared power laws of these Green functions. We discuss consequences 



0^ ' for the interpretation of recent results from lattice QCD. 
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I. INTRODUCTION 

In the past decade much progress has been made in the understanding of the low energy 
sector of QCD. This progress has been achieved both with continuum methods as well as 
lattice computations. In the continuum non-perturbative functional methods have been 
used: Dyson-Schwinger equations (DSEs) and functional renormalisation group equations 
(FRGs). Both frameworks are truly ab initio approaches in the sense that^he^can be 
derived rigorously from the full effective action of QCD, for reviews see (jj 
Although both frameworks constitute an infinite hierarchy of coupled equations, they a. 
for the extraction of scaling laws for Green functions in the deep infrared 
These scaling laws are related to important properties of the low energy limit of QCD, such 
as confinement and chiral symmetry breaking. 

A key building block relevant for the infrared behaviour of QCD are the ghost and gluon 
propagators. In Landau gauge the ghost dressing function gives access to the status of global 
gauge symmetry: an infrared diverging ghost unambiguously signals an unbroken symmetry 
corresponding to a well-defined global colour charge ^j]. This is an integral part of the 




low 
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Kugo-Ojima confinement scenario Q]. An infrared vanishing gluon propagator violates the 
Osterwalder-Schrader axiom of reflection positivity Q], and transverse gluons do not belong 
to the physical asymptotic state space of QCD. Finally, in Landau gauge one can construct 
a running coupling with a renormalisation group invariant combination of ghost and gluon 
dressing functions 

In Landau gauge QCD and in terms of correlation functions, the Kugo-Ojima confinement 
criterion is expressed as 



p 2 {A{p)A{-p)) 







p 2 {C{p)C{-p)) 



oo . 



(1) 



with the gauge field A and the ghost / anti-ghost fields C, C. An even stronger condition for 
the gluon propagator has been derived in a discretised version of Yang-Mills theory, where 
the infinite volume/continuum limit can be taken analytically: considerations on the impact 
re (first) Gribov horizon on dressing functions led to Zwanziger's horizon condition 



of t 

Q 



18] 



(A(p)A(-p)) 0, 



(2) 



which implies gluon confinement via positivity violation. 

The behaviour (P) and (j2J has first been seen in a DSE-study Q. This result has been 



confirmed and extended within further DSE-computations, e 



quantisation e. g. |9j , |22| as well as FRG-computations 
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24L l25j ]. for related work see 



In all these studies a non-renormalisation theorem for the ghost-gluon 



is used leading to 
p 2 (A(p)A(-p)} 



(P 2 ) 2K , p 2 (C(p)C(-p))^(p 2 )-\ (3) 

with k G [1/2, 1[. It has been argued in [l(3] that Q is the only consistent solution. Eq. (jSJ) 
has been extended to a self-consistent solution of the (untruncated) tower of DSEs in con- 
tinuum Yang-Mills theory 3|: for proper vertices Zq 
gluon legs the infrared asymptotics is given by 



as w ith n ghost, n anti-ghost and m 



7 (2n,m) f ^2 
J 0,as 
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where k is the exponent of the ghost dressing function (n = 1, m — 0) as defined in (jHJ). 
It is interesting to c omp are the continuum result (jSJ) with results from lattice QCD 
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the available finite volumes most lattice results for 
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the gluon propagator are compatible with JTJ and ©• Extrapolations towards the infinite 
volume limit, e.g. 3^|, seem to agree with an infrared finite propagator, i.e. 

p 2 (A(p)A(-p))~(p 2 ) 1 , (5) 

see however for an extrapolation leading to an infrared vanishing propagator. 

The situation is much less clear for the ghost dressing function. Whereas some simulations 
g ,ve an infrared diverging ghost Q Q, other anthers intent the, results as pointing 
towards an infrared finite ghost dressing function 38], i.e. 

(C(p)C{-p)) ~ (p 2 ) . (6) 

Clearly, (jHJ) and (0) together do not agree with the continuum result ©. Instead, they have 
been proposed as a second possible solution of the continuum DSEs [38]. 

In this work we shall show that the infrared asymptotics of Landau gauge Yang-Mills is 
uniquely fixed by DSE and FRG. We first discuss the relations between these two sets of 
equations in the next section. Then a general infrared analysis of the DSEs for the ghost and 
gluon propagators as well as for the ghost-gluon vertex is performed. In the following section 
we repeat this analysis within the FRG, and show that (0), © cannot survive the infinite 
volume /continuum limit. In section IS we show that self-consistency of DSEs and FRGs 
enforces the unique solution (j3J) for the infrared asymptotics of pure Yang-Mills theory. We 
briefly discuss the extension of the present analysis to QCD and the electro-weak sector of 
the standard model. In our concluding section we discuss the consequences of this result. 

II. FUNCTIONAL RELATIONS 

A quantum field theory or statistical theory can be defined uniquely in terms of its renor- 
malised correlation functions. They are generated by the effective action T, the generating 
functional of 1PI Green functions. For the present work we consider pure Yang-Mills theory 
with the classical gauge fixed action 

S c i = \j tiF 2 + ±J (d^f + J C-d.D^-C. (7) 

in the presence of an additional scale k, see (jj for a detailed discussion. The propagation is 
modified via fc-dependent terms 

AS k = ± [ A»R%Al+ [ C a R ab C b , (8) 
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where R a ^ v and R ab are k- dependent regulator functions. Within the standard choice k is 
an infrared cut-off scale, and the functions R cut-off the propagation for momenta smaller 
than k. Here we also consider more general R that only have support at about the mo- 
mentum scale k. Such regularisations allow for a scanning of the momentum behaviour of 
Green functions. The regularised effective action is expanded in gluonic and ghost vertex 
functions and reads schematically 

r * M = E -rii r f ,m) ° n ° n Am > ( 9 ) 

z — ' mini 2 

m,n 

in an expansion about vanishing fields = (A,C,C). In (jHJ) an integration over momenta 
and a summation over indices is understood. The effective action Tk satisfies functional 
relations such as the quantum equations of motion, the Dyson-Schwinger equations (DSEs); 
symmetry relations, the Ward or Slavnov- Taylor identities (STI); as well as RG or flow 
equations (FRGs). All these different equations relate to each other. Indeed, the Slavnov- 
Taylor identities are a projection of the quantum equations of motion, whereas flow equations 
can be read as differential DSEs, or DSEs as integrated flows. Written as a functional relation 
for the effective action T k and 8 pec,f yl n g to p Ur e YM theory, the DSE r ead S , ,, Q 

ST krn 6S c i r . , 



>opj , 



(10) 



where the operators op are defined as 



op (x) = / <?yG Hi [(j)]{x,y) — + <j>[x) , (11) 



and 



g -Mffs^L (12) 

is the full field dependent propagator for a propagation from 0i to 02, with = (A,C,C). 
The functional derivatives in (jlOj) act on the corresponding fields and generate one loop and 
two loop diagrams in full propagators. The functional DSE (|1(J|) relates 1PI vertices, the 
expansion coefficients of T^, to a set of one loop and two loop diagrams with full propagators 
and full vertices, but one classical vertex coming from the derivatives of S c \. We emphasise 
that the DSE ()10j) only implicitly depends on the regularisation via the definition of the 
propagator in (fE?|) . It has the diagrammatic representation 
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FIG. 1: Functional DSE for the effective action. Filled circles denote fully dressed field dependent 
propagators (|12j). Empty circles denote fully dressed field dependent vertices, dots denote field 
dependent bare vertices. 

Fig. ^ shows the structure of the functional DSE flT)|) . The rhs is given in powers of 
the field-dependent fully dressed propagator [</>], and its derivatives, as well as the field 
dependent bare vertices. The momentum scaling of Green functions is directly related to 
the scaling of these building blocks. 

The flow equation for the effective action reads 

dM<f>} = \Jd A P GM(P,P) dtR b ; u ( P ) - J d'p G ab [<P]{P,P) dtR ba (p) , (13) 

where t = Ink. The flow ()13j) relates the cut-off scale derivative of the effective action 
to one loop diagrams with fully dressed field-dependent propagators. We can contrast the 
diagrammatic representation of the DSE in Fig. ^ with that of (jlUjb 



d t r k [<j>] = 




FIG. 2: Functional flow for the effective action. Filled circles denote fully dressed field dependent 
propagators (fT2*)) . Crosses denote the regulator insertion dtR. 

Fig. El shows the structure of the functional flow (fT3|) . The rhs is given by the field- 
dependent fully dressed propagator G^cj)] and the regulator insertion d t R. Here, the mo- 
mentum scaling of Green functions solely depends on the scaling of G and d t R. If we choose 
the regulator function R such that it has the RG- and momentum scaling of the related two 
point function, the flow is RG-invariant and only depends on full vertices and propaga- 
tors, including the regulator term. The standard use of (|T3J) is to take a regulator function 
R(p 2 ) which tends towards a constant in the infrared and decays sufficiently fast in the 
ultraviolet, and hence implements an infrared cut-off. In the present context there is an- 
other interesting choice for R: let R only have support at momenta p about the scale k, and 
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R oc Tq ^ at momenta p 2 ~ k 2 . Then the regulator term does not change the theory, Y ^ ~ Tq, 
and ()13|) only entails the (change of the) momentum dependence of Green functions of Tq. 
This provides the resolution of the momentum dependence at p 2 of the full effective action 
To directly from the flow equation at k 2 = p 2 . 1 We shall detail this choice later. It is also 
worth noting that the relation between DSE and FRG is natural in a ATPI formulation 
which leads to a mixture of the scaling relations derived from and 1)13)) . 

In the present work we investigate the leading infrared behaviour of vertices and prop- 
agators constrained by consistency of (JTUJl and ([13)1 . The crucial ingredient in the related 
consistency equations is the fact, that the DSEs derived from (|T0|) also depend on bare or 
classical vertices whereas the flow equations derived from (|T3j) solely depends on full ver- 
tices. This allows us to extract non-trivial information of the theory from a finite set of 
vertex DSEs and FRGs that would require a whole infinite tower of either DSEs or FRGs if 
restricting the analysis to either of the functional equations. We analyse the leading infrared 
behaviour of (fTUj) and (fT^j) for momenta and cut-offs 

p\ k 2 « A 2 QCB . (14) 

To that end we introduce dressing functions Z^ n,m>l for one particle irreducible Green func- 
tions with n ghost, n anti-ghost and m gluon legs via 

if "' m) (pi, ..., P2n+m ) = Z^ m \ Vl , ...,p 2n+m )T^ m \ Pl , ...,p 2n+m ) . (15) 

The expansion coefficients r^ 2n,m ) Q f the effective action have been defined in ©. The 
7~(2n,m) denote the infrared leading tensor structure of the respective Green function, and 
carry their canonical momentum dimension. Then, following the IR analysis in |^ , the 
asymptotic vertex functions can be expanded about the leading asymptotics at vanishing 
cut-off k = 0: 

4 2 "' m) (Pl, :,P2n + m) * 4 2 a n s m) (Pl, J****) (l + ^^(Pl, -,P 2n+ni )) , (16) 

where pi = Pi/k, and the asymptotic infrared part Zq 2 ™^ only depends on ratios of mono- 
mials and possible logarithmic dependencies. Inserting the parameterisation ()16j) into the 



1 Such a single mode regulator cannot be used to solve the theory by successively integrating out degrees 
of freedom. However, it proves useful for studying fixed point solutions [4l| . 
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flow pUj) and solving for §Z^ 2n ' m ^ one can prove that 5Z( 2n ' m ) solely depends on pi. This 
suffices to fix the relations between the anomalous scalings of the vertex functions Z^™^ 
independent of the SZ^ 2n ' m \ 

For our analysis we only have to know the global scaling behaviour for the dressing 
functions Z^ n ' m \ that is, modulo logarithmic scaling, 

Z £' m) (M, Ap 2 n+m) = \ K2 ">™Z^ m \pi, ...,P2n+m) ■ (17) 

Specifically interesting for the Kugo-Ojima/Gribov-Zwanziger confinement scenario are the 
exponents /to,2 and ^2,0 °f the inverse gluon dressing function Z^ ' 2 ' , and inverse ghost dress- 
ing function Z^ 2 '°\ The horizon conditions (JTj),P| read 

«2,o > k ,2 < -1 • (18) 

We close this section with the remark on the interpretation of the scaling analysis derived 
from the combined functional relations (fTT1|) and (fT3|) . In principle, such an analysis produces 
the most singular scaling of all diagrams involved and is neither sensitive to cancellations 
between different diagrams nor to cancellations within a given diagram. However, to affect 
the infrared behaviour of the Green functions in a consistent way, such cancellations have to 
work in the whole tower of DSEs and FRGs and therefore can only be driven by symmetries. 
In the present case we consider such a possibility as highly unlikely. We will come back to 
these points at the end of section El 



III. INFRARED ANALYSIS OF GHOST AND GLUON DSES 

The Dyson-Schwinger equations for the ghost and gluon propagators are given diagram- 
matically in Fig. |3J The infrared behaviour of these propagators can be analysed as follows 
: We choose the external momentum scale p 2 , according to ()14J) . to be much 
smaller than any other scale, i.e. p 2 <C Aq CD , where A QCD ~ (9(200 Me V) is the non- 
perturbative scale of Yang-Mills theory generated via dimensional transmutation. The loop 
integrals on the right hand side of the DSEs are dominated by momentum configurations, 
where the internal loop momentum q is of the same order as the external momentum, i.e. 
p 2 ~ q 2 . The reason for this well known behaviour is the appearance of at least one propaga- 
tor in each loop with a denominator proportional to (p — q) 2 . Thus a self-consistent solution 
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of the equations with small external momentum can be obtained by replacing all dressing 
functions inside the loops with their infrared asymptotic behaviour. 






+ 




FIG. 3: Dyson-Schwinger equations for the gluon and ghost propagator. Filled circles denote 
dressed propagators and empty circles denote dressed vertex functions. 



We illustrate this analysis at the example of the ghost DSE. For the sake of comparison 
with the literature we switch to the standard DSE notation, where the non-perturbative 
dressing of the propagators is denoted by propagator dressing functions G{p 2 ) and Z{p 2 ): 

1 „, „ 1 



G( P 2 



4 2 '°V) 4°'V) 

and Zwanziger's horizon conditions (j!8j) read 

Z{p 2 ) 



Z(P 2 ) 



lim G(p ) = oo 



lim 

p 2 ^o p z 



The DSE for the ghost propagator reads is given by 

1 ~ ,„ f d A q G(q 2 )Z(l 2 ) 



Z 3 -g 2 N c 



P V{l)qZ^ l \p,q) 



(19) 



(20) 



(21) 



G(p 2 ) ° * V (2vr) 4 p 2 q 2 P 
with the momentum routing I = (q—p). The abbreviation pV(l)q denotes a contraction with 

(2 1) 

the transverse momentum tensor pV{l)q = p^P^^qu, and Zq ' (p,q) denotes the dressing 
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of the ghost-gluon vertex. The ghost renormalisation constant Z 3 absorbs all ultraviolet 
divergencies from the loop integral thus rendering the right hand side of the equation UV- 
finite. This can be made explicit within a momentum subtraction scheme. Here Z% is 
evaluated at a subtraction point p 2 = fi 2 , which we choose to be fi 2 = 0. One obtains 

*- 1 +r/N J«W rt(1 . ,22) 



3 G(0) * J (2vr) 4 4 g 4 
and subsequently 

gW) = W)~ 9 C J WVX vwv pV{l)qZ ° M 



3G(q 2 )Z(q 2 ^ 



Z^\0,q)Y (23) 



4 q A 

Now the integral is UV- finite and we replace the dressing functions in the loop by their 
infrared expansion in terms of the power laws 

Z{p 2 ) ~ (p 2 )- K0 ' 2 , G{p 2 ) ~ (p 2 )-^ , Z^\p, q) ~ {q 2 )^ , (24) 

The vertex function can be equally well represented by powers of (I 2 ) or, more realistically, 
by powers of (p 2 + q 2 + 1 2 ). The crucial point is, that after integration all powers of internal 
loop momenta will be transformed into powers of the only external scale p 2 for dimensional 
reasons. This can be seen easily for the expansion (|24|). which leads to integrals that can be 
performed employing 

{q 2 Y{l 2 f _ r(i + q)r(i + 6)r(-q-6) , b 

d<1 q 2 l 2 T(l-a)T(l-b)T(2 + a + b) {P) " (25) 

Plugging ()24j) into (|23p. performing the integration and matching with the left hand side, 
1/G(p 2 ) ~ (p 2 ) K2 '°, we obtain two self-consistent solutions: 

KP ' ' (26) 

(p 2 T ■ 

In the first case, the loop inte gral dominates the right hand side, and the constant 1/(7(0) is 
cancelled by other terms (see [lOj for a detailed discussion). In the second case this constant 
does not vanish and dominates the right hand side of the equation in the infrared. We thus 
end up with two possible conditions 

«2,o ~ < (27) 
0, 
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from the ghost-DSE. Either G(0) is divergent, in accordance with the horizon condition 
(pOjl. or it is finite as proposed in 13811 a nd violates (J20J). Strong arguments against the latter 
possibility have been discussed in |loi l^. where it has been concluded that k 2j o > 0. For 
the sake of the argument, however, we will not use this result here but proceed by exploring 
the consequences of both options. 

We would like to stress again that we could have obtained the solutions ()27|) without 
explicitly solving the loop integral: since the external momentum (p 2 ) is the only scale in 
our problem all powers of internal momenta in the loop have to translate into powers of 
external momentum after integration for dimensional reasons. Thus by simply counting the 
infrared exponents of all loop propagators and vertices we also arrive at ()27j) . 

We proceed by analysing the DSE for the gluon propagator. Schematically we can write 
this equation as 

^(jp,^ = %3 + n ta dpole(p 2 ) + n sunset (p 2 ) + H S quint(p 2 ) + ^-gluonloop{p 2 ) + ^-ghostloop{p 2 ) , (28) 

where the dressing loops appear in the same order as in Fig. |3J The static tadpole-term 
is absorbed in the process of renormalisation. We therefore have to analyse the infrared 
behaviour of the four remaining dressing loops. Counting IR-exponents in the loops we 
arrive at: 

IT„ t (p 2 ) ~ (p2)-3 K „, 2+ «o, 4 Ii squmt (p 2 ) ~ (^)-4«o, 2 +2ko,3 

Tlgluonloopip 2 ) ~ (p 2 )- 2 ^+^ U ghostloop ( P 2 ) ~ (/)- 2 ^0+« 2 ,l ( 29 ) 

The infrared leading term from the right hand side has to match the left hand side 1/Z(p 2 ) ~ 
(p 2 ) K0 > 2 of the DSE. We thus obtain the expression 

k ,2 = min(o, -3k , 2 + k ,4 , -4k 0i2 + 2k , 3 , -2k ,2 + ^o,3 , -2k 2 ,o + «a,i) • (30) 
and subsequently 

k 0j2 = min(0, ^k ,4 , g«o,3 , ~/c 0) 3, -2k 2 ,o + «2,i) , (31) 

as our final condition for the gluon exponent kq^ from the gluon-DSE. In general we expect 
«o,2 < according to the Kugo-Ojima and horizon conditions (JTJ), (f2"U|) and lattice QCD 
(see e.g. ^ * s interesting to note that for negative k ,2 the contribution from the 

gluon loop, 1/3^0,3, is never the leading one on the right hand side of (|5D|). since it is always 
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FIG. 4: Dyson- Schwinger equation for the ghost-gluon vertex, derived via Eq. (|lt)|). All internal 
propagators are taken to be fully dressed. The ellipsis denotes two-loop diagrams, which are not 
needed for our analysis. 

dominated by the contribution from the squint diagram, 2/ 5^0^ Thus any truncation of the 
gluon-DSE that assumes a leading gluon loop (see e.g. |44J, l45l l46j] ) is missing the dominant 
contribution in the infrared. 

A further crucial ingredient is the DSE for the ghost-gluon vertex. One version is derived 
from the DSE for 5T/SA, and is given diagrammatically in Fig. 0] Similarly to the ghost 
and gluon propagator DSE we arrive at 

k 2 ,i < min(2K 2i0 + k 0i2 , 2k 0i2 + k 2 ,o , ^2,2 - 2« ,2 , ^4,o - 2« 2> oJ • (32) 

The inequality takes into account that the exponents of the two loop diagrams in the vertex- 
DSE may even be smaller than those of the one-loop diagrams considered in (J32j) . 2 From 
(J22J) we conclude that 

«2,i ^ 2k 2j0 + k ,2 • (33) 

Eq. (JHHj) together with the FRG-relation derived in the next section suffices to uniquely fix 
the relations between all Km,m in a closed form. 

Based on the conditions (J27j) and (|31|) and the exact equality in (|33|) an infrared analysis 
of the DSEs for the three-gluon vertex and the four-gluon vertex has been performed in 



12j. These results have been generalised to any Green function with n external ghost, n 



anti-ghost and m gluon legs: 

CV)-^". (34) 

with k = /t 2) o > 0. This expression solves (|27j) . (|ST|l and any other condition from the 

higher DSEs. In addition it solves the Slavnov- Taylor identities. Important aspects of this 



2 K2,i can be also determined from the DSE for 5T/5C or 5T/8C, see Fig.^ However, the present analysis 
then turns out to be more complicated even though two loop terms are absent. 
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solution are discussed in detail in J3J. Two of the characteristic properties of (J54*j) are: 0) 
contributions from ghost-loops always dominate the DSEs and (ii) it leads to IR-fixed points 
in the running couplings from the ghost-gluon (gh), three-gluon (3g) and four-gluon vertex 
(4g). These couplings are defined via 

a*V) = ^[Z^ l \p 2 )fG 2 {p 2 )Z{p 2 ), (35a) 

« 3 V) = f^'V^V), (35b) 

a 4 V) = V)] ^V) , (35c) 

where Z^ 1 = Zq 0,2 -* and G _1 = Zq 2 '°\ The vertex dressing functions Zq 2 ' 1 ' ) (p 1 ,P2,P3), 
Zo°' 3 \pi iP2,P3) and Zq°' 4 ^(pi,P2,P3) are evaluated at the symmetric momentum point 
Pi = v\ = V% = V 2 1 which make them functions of p 2 only. 

From the tower of DSEs alone it is difficult to prove, that (}3"4*)) is unique. One way to 
search for a second possible solution would be to assume K2,o = and /to,2 = — 1 from the 
start, corresponding to the behaviour © and (0), as proposed in j^. From Eqs. (J27)l and 
(f3T) one obtains consistency provided one of the three vertices is strongly divergent: 

k 2 ,i = -1, or k ,3 = -5/2, or k oa = -4 . (36) 

In the next section we will show that all options lead to inconsistencies in the 
functional renormalisation group equations. 3 As discussed in section |HJ any solution of 
the tower of DSEs necessarily has to solve the tower of FRGs as well. This provides tight 
constraints on possible solutions, which are in fact sufficient to prove the uniqueness of ()34j). 
as we shall see. 

IV. INFRARED ANALYSIS OF GHOST AND GLUON FLOWS 

Now we repeat the infrared analysis within the FRG framework. We restrict ourselves to 
regulator functions of the form 

R k (p 2 ) = T {2 \p 2 )r(p 2 /k 2 ), (37) 



3 Note that the second option together with l|35b|) leads to a strongly divergent running coupling, a 3g ~ 1/p 2 , 
which appears to the lattice results of 01- 
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where is the corresponding two point function T^' ^ for the ghost, and T^' 2 ^ for the gluon. 
Regulator functions (|3~7f) guarantee the persistence of the standard RG-scalings in the pres- 
ence of an IR cut-off, and are best-suited for the present studies. Within the parametrisation 
(fr7)|) and as a consequence of (fHTj) propagators G{p 2 ) take the asymptotic form 

G{p2) = ifV) i + W) ~ rg s (^)i + ^( 2 )(p)i + r(^) • ( 38 ) 

Eq. ()38)1 can be solely written in terms of p and /c-dependences. Then it reads 

/ r " 1 1 (ffl) 

r^(p2)l + ^)(p)l+r(p 2 )' 1 ) 

where k is either k ,2 (gluon) or k 2 ,o (ghost). The same rescaling can be done with all vertex 
functions: 

rf 1 ..., P2n+m ) ~ at**- r£' m) (pi, -,P2n + m)(i + ^ (2n - m) (Pi, Ate+m)) • (40) 

Another option for the infrared analysis is the choice of a regulator function Rk that only 
has support for momenta at about k 2 : 

R k (p 2 )=T {2 \p 2 )6 e (p 2 -k 2 ) (41) 

where 5 e (x) is proportional to a smeared out 5-function at x — 0, the <5Z( 2ra ' m ) only have 
support at momenta p 2 ~ A; 2 . With regulators (|41jl . the momentum dependence of r^ 2n,m ) 
agrees with that of r^ 2n '\ Only the strength of the two-point function rjj, ~ in the 
momentum window p 2 rs k 2 is changed. In particular, the infrared power laws at k ^ 
agree with those at k — 0. Therefore we can directly read-off the momentum-dependence at 
the momentum scale p 2 = k 2 . 

In turn, the general analysis with (|4()jl provides additional information on the infrared cut- 
off flow. Within the parameterisation (f4TJJ) integrated asymptotic flows for general vertices 
read 

sz^fa, jw») - f ^ ( ^y^rf -A (& ...,p' 2n+m ) , (42) 

^° V 0,as / 

with possible sub-leading terms. Eq. (|42J1 defines consistently renormalised finite DSEs 
In contrast to the DSEs (fTUj) it solely depends on full vertices but also involves an integration 
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over the cut-off scale k. The term dtT <y2n,m ^ on the rhs of (|42jl is derived by taking 2n,m- 
derivatives of the rhs of the flow (|13|). leading to a sum of one loop diagrams with dressed 
vertices and dressed propagators. For vertices with K2n,m < we have 

<5Z (2n ' m) (0,...,0) = -1. (43) 

Eq. ()43j) simply entails that an infrared cut-off is present and the divergent infrared behaviour 
for k = is suppressed. Trivially the infrared limit (|43J) only depends on p. From (}4*3*j) we 
derive a relation between the involved re^- with i,j < n + 1, m + 2 in the flow of r^ 2n,m ) 
within an iteration about <5Z( 2n ' m ) = on the rhs of (}4~2*j) . The analysis for K 2 n,m > within 
the integrated flow ()42j) is a bit more involved in one to one correspondence to possible 
difficulties with bare terms in the DSEs. However, as outlined above we can also directly 
resolve the momentum behaviour from (fT3*|) with regulators (jH$ , where these problems are 
absent. In fact this eliminates the possibility of solely dominating bare terms. 

We proceed with the analysis of the propagator FRGs. They can be derived from Fig. El 
and are shown diagrammatically in Fig. 

We exemplify the analysis at the gluon propagator with /to,2 < 0. As SZ^ ' 2 ^ has 
to approach —1, the momentum scaling of the k- integral has to precisely cancel that of 
l/T ( Q> 2 \p -»• 0). On the rhs of flU we can iterate the full Y^ n ' m ^ about those at k = 0, 
rg^ ,m ', Consequently we can simply sum over the H2n,m to identify the leading p-behaviour. 
With the regulator (}4*T]) this follows directly, as within this choice we have Yf n ' m ^ oc 
in the flow. Then the relations between the K2n,m follow from simple counting of powers of 
momenta. From the flow equation for the gluon propagator we derive the relation 

K 0j2 = min^2K 2 ,i - 2k 2 , , 2k , 3 - 2k , 2 , «o,4 - «o,2, «2,2 - «2,o) , (44) 

which can be solved for « 0) 2, 

k ,2 = min^2/t 2 ,i - 2k 2 ,o , |«o,3 , |«o,4, ^2,2 - «2,qJ • (45) 

Eq. (j43j) is already sufficient to rule out all three options in (|3Tj) . Indeed, if we insert (|3T)j) 
into (J43j) we arrive at 



-2 



«0,2 < < 



I (46) 
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FIG. 5: Functional renormalisation group equations for the gluon and ghost propagator. Filled 
circles denote dressed propagators and empty circles denote dressed vertex functions. Crosses 
indicate insertions of the infrared cutoff function. 

for the three options. However, goes with n 0>2 = — 1. The behaviour (jSJ) and (jBJ), 
proposed in |38j, is therefore ruled out. 

Now we use the combined DSE-FRG analysis to uniquely determine the re's without any 
further input. We shall see that a self-consistent solution leads to 



re 2 ,i = 



(47) 



To that end we also discuss the derivation of the re's for the ghost-propagator and the 
ghost-gluon vertex. For general regulators the infrared analysis for the ghost propagator is 
intricate and we defer the reader to [ill . I25I ] . With the choice (141 J) the result follows directly 
from the flow of the ghost propagator. Analogously to (J43j) we get from Fig. 



1 1 

K 2fi = min( K 2 ,l - -K ,2 , «2,2 - K 0i 2 , -«4,0 



(48) 



The FRG relations (J43j) . (}4*%|) as well as the DSE relations (f2Tj) . (p?Tj) for ^2,0 and ^0,2 are not 

closed as they also depend on vertex kappa's. The ghost-gluon vertex comprises the crucial 
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information. It is protected by non-renormalisation which turns out to be powerful enough 
to fix the whole system completely. Its flow is given by all one loop diagrams with regulator 
insertions and full vertices (up to 5 point vertices) with one external gluon line and one 
ghost and anti-ghost line. It reads schematically 



kdk 



.6 




(a) 



XX 



(d) 



+ 



(b) 




(c) 





(e) 



(/) 



is) 



FIG. 6: Functional renormalisation group equations for the ghost gluon vertex. All propagators and 
vertices are fully dressed. Only one possible insertion of the infrared cutoff function per diagram 
is shown. 



The infrared analysis of Fig. El with (|42|) . or alternatively employing (|41jl. leads to 
« 2) i = min^2/c 0)2 + «2,o - ^0,3 , ^0,2 + «>2,o , ^0,3 + ^2,2 - 2k ,2 , «4,i - k 2 ,o , k 2 ,3 - ^0,2) ,(49) 

from the diagrams (a), (6), (c), (/), (g) respectively. The diagrams (d),(e) involve exactly one 
ghost gluon vertex and the related anomalous scaling cancels on both sides of Fig. El The 
surviving k's cannot be negative, as this spoils self-consistency of the integrated flow. Hence 
(d),(e) lead to the constraints 



«4 i0 — 2k 2 ,o > 

^2,2 - K-0,2 - «2,0 > 



(50a) 
(50b) 



Additionally the non-renormalisation of the ghost-gluon vertex [20] constrains 



«2,1 < 



(51) 



as at least one tensor structure of the ghost-gluon vertex has a finite dressing. From the 
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second term on the rhs of (|49|) we extract 

«2,1 < |«0,2 + «2,0 • (52) 

We insert (JH2|) in the first term of (f4*4"j) and arrive at k,q,2 < ^0,2- Consequently the bound in 
f!52j) has to be saturated and 

«2,1 = 5«0,2 + ^2,0 • (53) 

The DSE-analysis leads to the constraint /c 2) i < ^0,2 + 2k 2 ,o> ff55f) . With (J53|l this turns into 
1^2,1 < 2«2,i- As ^2,1 < 0, (JSH), we arrive at the unique solution 

« 2l i = , (54) 

accompanied by the relation 

k 0)2 = -2k 2 , , (55) 
for the dressing of ghost and gluon propagators in agreement with jsj]. 



V. UNIQUENESS OF INFRARED ASYMPTOTICS 

The analysis of the last two sections for the propagators and the ghost-gluon vertex can 
be extended to all K n>rn . We first derive the relations for the purely gluonic three and four 
point functions. The flow of three and four gluon vertices is given by all one loop diagrams 
with regulator insertions and full vertices (up to 5 and 6 point vertices respectively). For 
the three gluon vertex this reads schematically 

Using Fig. [7| we arrive at 

k 0j3 = min^|K ,2 , 3« 2 ,i - 3k 2 , , ^2,2 + «2,i - 2k 2 ,o , ^2,3 - ^2,0 , ^0,5 - «o,s sj , (56) 
from the diagrams (a), (b), (d), (e), (/). The diagram (c) leads to the constraint 

«o,4 - 2k ,2 > , (57) 

which already restricts the singular behaviour of the four gluon vertex. We also remark that 
(JSnj) puts a simple upper bound on k ,3, namely 

K 0,3 < § K 0,2 • (58) 
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FIG. 7: Functional renormalisation group equations for the three gluon vertex. All internal prop- 
agators are taken to be fully dressed. Only one possible insertion of the infrared cutoff function 
per diagram is shown. 

This bound is the natural scaling of the vertex in the presence of a fixed point for the 
coupling constant a s . Note however that a priori not all couplings as defined in (JHo]) run to 
a fixed point. 

The same analysis can be done for the four gluon vertex. Its flow reads schematically 



kdi 





(6) 





(rf) 



i 




FIG. 8: Functional renormalisation group equations for the four gluon vertex. All internal propa- 
gators are taken to be fully dressed. Only one possible insertion of the infrared cutoff function per 
diagram is shown. 



Similarly as for the three gluon vertex we derive from Fig. |S] the anomalous scaling of the 
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four gluon vertex 

k ,4 = min^4/t ,3 - 4k ,2 , 4/c 2) i - 4k 2 ,o , «o,5 + «o,3 - 2k , 2 , ^2,3 + «2,i - 2k 2 ,o , 

2k ,2 , 2/« 2 ,2 - 2/t 2 ,0 , ^0,6 - «0,2 , «2,4 ~ «2,0 , «2,2 + 2^2,1 ~ 3/€ 2 ,0^) , (59) 

from the diagrams (a), (h), (j) respectively. The diagram (i) leads to the constraint 

2k ,3 - 3k ,2 > . (60) 
The first term in the second line of (p)9"j) puts a bound on k ,4; 

^0,4 < 2k ,2 • (61) 
Together with the constraint (fBTj) this gives the unique solution 

k oa = 2k 0j2 , (62) 

and 

«0,3 = §«o,2 • (63) 

Note that the scaling laws ()62)) . (j63)) are already generated by the diagrams involving ghosts. 
Indeed, this is valid for all proper vertices and leads to the unique solution 

«2n,m = {n — m)K, with k — k 2 ,o, (64) 

which we now prove: first we observe that with (|64|) all diagrams in the FRG have the same 
leading infrared asymptotics. Let us assume for a moment that (|64j) is not true for all proper 
vertices. Then at least one vertex has 

«2n ,m < { n - m )K . (65) 

The vertex r( 2n °> m °) occurs in diagrams of FRG for lower vertex functions Y <y2n,m " 1 with 
no — n = 1 or mo — m = 1. Necessarily also these vertices satisfy ()65p and disagree with 
(164)) . Within an iteration this enforces that all vertices with n < Uq and m < m satisfy 
flfiljj) . in particular /t 4 , k 3 , k ,2, «2,i> K 2,o- This contradicts the uniqueness of the results (jS3), 
flo3|) .()6^ |) .()6^ |) derived above, and hence proves (jMl) . 



For the special case of the propagators this relation has been already derived in 



with the help of additional physical constraints. In |l2j the self-consistency of (I64j) for the 
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whole tower of vertex DSEs has been shown. The dominance of ghost loops in the tower of 
DSEs is equivalent to the dominance of the Faddeev-Popov determinant over the Yang-Mills 

n 

action, as proposed in 22]. In the present work we were able to extend these results to a 
proof of uniqueness based on a self-consistency analysis of the quantum equations of motion. 

The above results hinges on a key structure valid for general theories in the presence of 
a single dynamical scale, and follows already from the structure of the functional DSE (|10jh 
Fig. ^ and FRG (jl3j) . Fig. O any vertex DSE comprises a sum of diagrams proportional to 
a subset of the bare vertices of the theory at hand. Consistency with the FRG, which only 
depends on dressed vertices, requires that in all vertex DSEs at least one of these vertices, 
if dressed, has ^vertex = 0. In Landau gauge Yang-Mills this is the ghost-gluon vertex. In 
general the above criterion leads to more than one vertex with fiwtex = 0. 

As an example for this general pattern we extend the pure gauge theory analysis to 
a YM-Higgs theory. The functional DSE and FRG can be read-off from (|lUj h (|13j ) with 
= (A, C, C, h) and an additional Higgs action Shi ggs = \j(D ll h) 2 + V[h]. Under the 
assumption of a single, dynamical mass scale we deduce a unique relation for a general 
vertex function r( 2ra,m,2/ ) with n ghost, n anti-ghost, m gluon, h Higgs lines, 

«2n,m,/i = {jl - TTl)K with K = K 2 ,0,0 ■ (66) 

In particular it follows that the Higgs propagator has a constant dressing: = ^0,0,2 = 0, 
as well as the 4 -coupling «o,o,4 — required by the presence of vertices with constant 
dressings in all DSEs. Note that K^n.m,/ — includes logarithmic scaling. Eq. (Jol))) is only 
valid in the symmetric phase. In the spontaneously broken phase we expect a massive gauge 
field propagator, ^0,2,0 = ~ 1, and massive Higgs propagators, = —1. Furthermore, as a 
positive k > for the ghost signals an unbroken (colour) symmetry, we conclude k < in 
agreement with the converse of the Higgs theorem Due to the additional mass scale 
some if not all other vertices may scale canonically, i.e. K2 n ,m,h — 0. The present infrared 
analysis then shows consistently H2n,m,h > 0, no singular scaling occurs. This is in marked 
contrast to massless Yang-Mills theory. 

In the case of full QCD, including dynamical quarks, the present infrared analysis has 
interesting consequences which shall be published elsewhere. Finally we discuss the caveat 
mentioned at the end of section |Hj From the above analysis it is clear non-trivial cancel- 
lations always have to occur in an infinite sub-set of diagrams. It is hard to see which 
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symmetry should be responsible for such a behaviour, as constraints from gauge symmetry, 
i. e. STIs, are respected by the solution . 

VI. CONCLUSIONS 

In this work we used Dyson-Schwinger equations (DSEs) and functional renormalisation 
group equations (FRGs) to analyse the infrared behaviour of proper vertices of SU(N C ) 
Yang-Mills theory. We have shown that the structure of these functional relations is suf- 
ficiently different to generate tight constraints for infrared anomalous dimensions of these 
vertices. The caveats of this construction have been discussed at the end of sections |H] 
and The constraints are powerful enough to enforce a unique solution for the 

infrared behaviour of proper vertices in the presence of only one external scale. Thus the 
Kugo-Ojima criterion is satisfied ensuring a well-defined global colour charge. A further 
consequence is the fixed point behaviour of the running coupling in the infrared, since this 
behaviour is implied by the solution (@J) via the non-perturbative definitions in Eqs. (jSHJ)- 
In turn the proposal (J5J), (0) in j3] is excluded. 

We emphasise that both, the similarities as well as the differences of DSEs and FRG were 
crucial for our results. This structure is certainly useful beyond the present investigation, 
for example if devising truncation schemes. Implicitly this was already used for enhancing 



the respective reliabi 



ity: the coinciding results for the infrared asymptotics from DSE 




0, |10|, |l2j and FRG l23l . l25j ] are non-trivial as the functional equations are sufficiently 
different. 

The present consistency analysis not only uniquely fixes the infrared asymptotics but 
also excludes certain truncation schemes of DSEs and FRGs: e.g. we have shown that 
truncation schemes of the gluon-DSE in Landau gauge Yang-Mills relying on an infrared 
leading behaviour of the gluonic one-loop diagram, as e.g. assumed in |3] miss the 
leading infrared behaviour. In addition, truncation schemes that assume all terms in the 
g ,uon-DSE to be e^aUy leadi„ g are excluded as well. 

It would be desirable to reproduce (@J) from lattice QCD. To this end one has to ad- 
dress some caveats in comparing infrared results from the lattice with those of a continuum 
approach. Lattice simulations are necessarily performed at finite volume and finite lattice 
spacing and one has to carefully perform both an infinite volume and continuum limit extrap- 
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olation. These procedures are currently under debate [37|, |39j, |48|, |49j, |51( . It is interesting, 
however, that the procedure of [stJ gives Kq^ ~ —1-04 in agreement with (^J). Unfortunately 
direct lattice calculations in the infrared scaling region p < 100 MeV are extremely expen- 
sive in terms of CPU-time and have not yet been performed in four dimensional Yang-Mills 
theory. This is different in three dimensions, where lattice results are in good agreement 
with the corresponding power-law analysis in the continuum jsoj l . 

Furthermore, gauge fixing is implemented differently In the continuum theory one uses 



either the Faddeev- Popov method 



22| ) , whereas on the lattice 



52j or stochastic gauge fixing 
a gauge fixing functional is extremised. Due to the presence of Gribov copies this might affect 
the infrared behaviour of Green functions. These effects are currently under investigation in 
the continue □ and on the lattice QQSHHQ. The effects 8 ee,n to be ,nnch 
stronger for the ghost than for the gluon propagator. This nicely corresponds to the fact 
that lattice and continuum solutions agree much better for the gluon than for the ghost. 

The present analysis can be extended to full QCD, and reveals an interesting structure. 
Related results will be published elsewhere. 
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